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Motivated by the observation of light surface states in SmBe, we examine the effects of surface 
Kondo breakdown in topological Kondo insulators. We present both numerical and analytic results 
which show that the decoupling of the localized moments at the surface disturbs the compensation 
between light and heavy electrons and dopes the Dirac cone. Dispersion of these uncompensated 
surface states are dominated by inter-site hopping, which leads to a much lighter quasiparticles. 
These surface states are also highly durable against the effects of surface magnetism and decreasing 
thickness of the sample. 


Kondo insulators are a class of strongly correlated elec¬ 
tron material in which the screening of local moments 
by conduction electrons gives rise to an insulating gap 
at low temperatures [Ilia. The first Kondo insulator, 
SmBe, discovered more than 40 years ago [3], has at¬ 
tracted renewed interest due to its unusual surface trans¬ 
port properties: while its insulating gap develops around 
Tk — 50K, the resistivity saturates below a few Kelvin[4]. 
Although this excess conductivity was originally ascribed 
to mid-gap impurity states [5], renewed interest derives 
from the proposal that SmBe is a topological Kondo in¬ 
sulator, developing protected conducting surface states 
at low temperatures [BHS] . Experiments [TUHTS] have 
since confirmed that the plateau conductivity derives 
from surface states, and these states have been resolved 
by angle-resolved photoemission spectroscopy (ARPES) 
[HHn]. The most recent spin-ARPES experiments have 
also resolved the helicoidal spin polarization of the sur¬ 
face quasiparticles expected from topologically protected 
Dirac cones [18]. 

Yet, despite this success, certain aspects of these ma¬ 
terials remain unexplained. A particularly notable prob¬ 
lem, is that both quantum oscillation [19] and ARPES [TSE 
[T8] studies show that the surface quasiparticles are light, 
with Eermi velocities Vg, ranging from 22QmeVA[T4] to 
300 meVA [15]; yet current theories [BHS] predict heavy 
Dirac quasiparticles with velocities Vg 30 — 50 meVA, 
an order of magnitude slower. 

Here we propose a resolution to this problem, by tak¬ 
ing into account of the breakdown of the Kondo effect 
at the surface. The essence of our theory is based on 
the observation that the reduced co-ordination of the 
Sm^+ ions at the surface causes a marked reduction in 
the surface Kondo temperature ^ Tx/lO, so that the 
screening of local moments at the surface is either sup¬ 
pressed to much lower temperatures, or fails completely 
due to an intervention of surface magnetic order. “Kondo 
breakdown” [20142^ liberates unquenched moments at the 
surface, and has the effect of shifting the mixed valence 
of the surface Sm ions towards the higher entropy, 4/^ 
(3+) configuration of the unquenched moments, as ob¬ 
served in X-ray absorption spectroscopy [25]. Most im- 
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FIG. 1: (a) Schematic showing the the non-local Kondo in¬ 
teraction at the surface of topological Kondo insulator. The 
opposite parities of the / and d bands cause the hybridizaBon 
to acquire a p-wave, spin-dependent form factor 4>(R) oc iR-a 
(see eq. |^, which vanishes onsite, (b) and (c) respectively 
contrast the surface states of a topological Kondo insulator, 
before and after surface Kondo breakdown, showing how the 
Dirac point drops into the valence band, giving rise to lighter, 
faster surface states. 


portantly, Kondo breakdown will liberate a large number 
of carriers previously bound inside Kondo singlets at the 
surface (see Eigl(b) and (c)): this process localizes 4/- 
holes, driving up the electron count in the surface states 
to form large spin-polarized Eermi surfaces determined 
by the Luttinger sum rule 


(27r)2 


An/ 


( 1 ) 


where Mfs' is the total Eermi surface area of the singly- 
degenerate surface states, Auf is the change the / va¬ 
lence, i.e. for Sm 4/^-^+ ^ 4/^+, An/ = 0.4, in units 
where the lattice constant a = 1. A detailed analysis 
presented later, shows that the dispersion of these highly 
doped surface states are dominated by a quadratic inter 
site hopping term 


E{kx, ky) = Eo + ^jTktfr{5)k^ + tfp{S)k]_, (2) 


where t/ is the effective / hopping, while r{S) and p{S) 
are functions that depend on the surface scattering phase 
shift 6, as defined later. The quadratic term, absent in a 
conventional topological surface state, appears when the 
perfect compensation between heavy and light electrons 
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is disrupted at the surface by Kondo breakdown: this 
term increases the velocity of the surface states by about 
an order of magnitude, providing a natural explanation to 
the light surface states observed in experiments. More¬ 
over, the heavy doping of the surface states shifts the 
Dirac point into the valence band, leading to an unusual 
protection of surface states against decreasing thickness 
and surface magnetism. We will show that the interac¬ 


tion of the local moments and the light surface states is 
described by a new kind of “chiral Kondo lattice” with 
the potential for a rich phase diagram of competing in¬ 
teracting surface states. 

Model: We introduce a simplified lattice model for topo¬ 
logical Kondo insulators, described by the periodic An¬ 
derson model: 


H = 


E (4s. 4) 
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where the and respectively create conduction and 
f electrons at site j, with corresponding hopping matrix 
elements —tij and i{j . U is the onsite Coulomb repulsion 
of f-electrons. V and Cf are the bare hybridization and 
/ level position. The form-factor 


^(^)/37 = 


n.n. 

0 otherwise 


( 4 ) 


describes the spin-orbit coupled hybridization between 
neighboring / and d electrons. This simplified form fac¬ 
tor captures the important vectorial spatial structure of 
hybridization between f and d states, whose orbital angu¬ 
lar momenta differ by one unit. The hybridization is odd- 
parity and vanishes onsite because the heavy and light 
orbitals have opposite spatial parity. Fig. 1(a) shows a 
schematic of our model. Note that it explicitly includes 
an / dispersion [2], which is required for a fully gapped 
spectrum. Microscopically these terms are expected to 
result from the indirect hopping of f-electrons via the p- 
orbitals of the Boron ions. 

To describe the large U limit we adopt a slave-boson 
mean-field theory and carry out a saddle point approxi¬ 
mation in the bulk, giving rise to the mean-field Hamil¬ 
tonian [26] . 
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Here, V = Vb is the renormalized hybridization, where b 
is the slave boson projection amplitude. The f-hopping 
now becomes = b‘^i-^. We consider a general near¬ 
est neighbor and next nearest neighbor dispersion, Ck = 
—2t cos ki — 4t' ^i^j cos ki cos kj while e/k = 
where, for simplicity, we have taken the ratio between 
the /— and d— electron hoppings to be a single fixed 
constant a = tf /t = for both nearest neighbor and 
second-neartest neighbor hoppings. The quantity A is the 


constraint field that imposes the mean-field constraint 
Q = Uf +6^, where Q is the local conserved charge as¬ 
sociated with the slave boson treatment of the infinite 
U limit, here taken to be Q = 1; A4 is the total num¬ 
ber of sites. In momentum space the form factor of the 
hybridization now takes the form 4>(k) = 5k • a, where 
5k = (sin kx , sin ky , sin kz ), which reduces to k • <j at small 
k. This odd-parity hybridization is reminiscent of the gap 
function in topological superfluid ^He-3B. This model can 
be regarded as an adiabatic continuation from small to 
large U at infinite spin-orbit coupling. Band crossing 
between the odd and even parity bands generates the 
topological band structure [27l[28] with protected surface 
states. 

In order to obtain the surface state spectrum, we adopt 
Volovik’s approach, mapping the reflection at the bound¬ 
ary onto transmission through an interface where the hy¬ 
bridization changes sign[T]. This method enables us to 
solve for the surface states using a linearized Hamilto¬ 
nian. For pedagogical purposes, we take t' = = 0. 

We can treat the surface eigenstate at = 0 as a one¬ 
dimensional problem, and since the main surface scatter¬ 
ing effect takes place close to the Fermi surface, we can 
linearize the dispersion([^ normal to the surface and ob¬ 
tain the surface eigenstates at kx^y = 0. The two interface 
eigenstates at kx^y = 0 are given by 
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where kz = X/[t{a + 1)] and k = V/{^/at). r and cf 
denote the orbital and the spin component of the wave 
function. The corresponding energy of a state at the 
Dirac point is Eq = A/((a + 1). The transverse dis¬ 
persion that develops at finite kx^y is then treated us¬ 
ing first order degenerate perturbation theory, projecting 
the full Hamiltonian Hmf onto the surface bound-states. 
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Width of the slab 


FIG. 2: Slab calculation for topological Kondo insulator: (a) for the uniform solution, (b) surface Kondo breakdown: surface 
local moments decouples from the rest of the system. Notice the large Fermi surface with high velocity quasiparticles in the 
case of the surface Kondo breakdown, (c) Inhomogeneous mean field theory calculations for the slave boson amplitude, 6 as a 
function of the thickness of the slab, b is reduced at the surface and for 'J. Tk, the local moments completely decouple. 


Heff = {i^±\HMF\i^±)‘ This leads to the dispersion 
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where kj_ = y+ k‘^. Since a <C 1, the surface states 
are composed of heavy quasiparticles. Note that at this 
stage, the quadratic k‘^ term is absent. 

Surface Kondo breakdown: The reduced co-ordination 
number of the f-electrons, plus the internal gapping of 
the bulk states lowers the effective Kondo coupling con¬ 
stant. If we assume crudely that the Kondo coupling 
constant is reduced by a factor of 5/6th, the reduced 
Kondo temperature is then 

T, = De-^.^Ti = = T, 

(8) 

For example, if we take the bulk Kondo temperature 
Tk ^ 50iF with half band-width D ^ 2eV ^ 2 x lO^iF, 
this gives a surface Kondo temperature of ^ 15iF. 
In practice, the gapping of the bulk states makes this an 
upper bound. This effect has been recently explored in 
1D[301 [3l]. However, unlike one dimension, where the 
spin-decoupled phase is no longer topological [30], in two 
and three dimensions, the topological insulator is pro¬ 
tected by time reversal symmetry and the spin-decoupled 
state remains topological in the bulk. With these consid¬ 
erations in mind, we model the surface Kondo breakdown 
as the suppression of the slave boson amplitude b to zero 
on the surface layer of the TKI. 

Surface Kondo breakdown localizes mobile /- 
quasiparticles at the surface, removing them from the 
Fermi sea. In SmBe, this is equivalent to a shift in the 
effective f-valence from 4/^-^ ^ 4/^ at the surface. This 
is a reduction in the number of f-holes, that corresponds 
to an increase An/ = 0.4 in the density of electrons 
on the surface. From Luttinger’s sum rule (eq. [^, the 
area enclosed by the Fermi surface is equal to the total 
electron density. The Fermi surface area observed in 
ARPES is 0.3-0.35 [T3UT8] . in agreement with Luttinger 
sum rule. 


Next we examine the effect of Kondo breakdown on 
the dispersion of the surface states. The effect of Kondo 
breakdown can be incorporated by introducing a scatter¬ 
ing phase shift 5 = kzU into the / component of the sur¬ 
face state wave function, which ensures that surface node 
of the / part of the wave function is shifted one lattice 
unit in towards the bulk, depleting the surface of mobile f- 
electrons (see supplementary information for details [32] E 
In Volovik’s approach, 6 becomes a transmission phase 
shift. The unperturbed energy of the surface bound-state 
is now modified Eq = (A + 2Vy/a sin 6)/{1 + a) and the 
dispersion becomes 

2k5 

E'{k,, ky) =E'^ + VV^r{5)k^ + at ^ ^ kl (9) 

ae +1 

Note the appearance of a k\ term, driven by an in¬ 
crease the conduction electron content of the surface 
states. The subleading linear term is also modified, con¬ 
taining the renormalization r{5) of order one: r{5) = 
2(cos(5-h ^ sin(5)/(l+ aexp(^)). Substituting at = tf 
and V ^ ^/TkI in eq. [pleads to eq.|^ This analytic result 
provides a natural explanation to the light quasiparticles 
in SmBe. For a parameter choice of t = 500 meV, tf = 
5 meV, consistent with band structure calculations [9], 
and the value of V = 100 meV which leads to the exper¬ 
imentally observed gap A 20 meV we obtain surface 
state velocity Vss ^ 58 meVA without surface Kondo 
breakdown and Vss ^ 340 meVA with surface Kondo 
breakdown in agreement with the Vgs ^ 300 meVA mea¬ 
sured by ARPES experiments [15]. Here the quoted sur¬ 
face velocities were calculated at the mid-gap energy of 
the Dirac point in the absence of Kondo breakdown. We 
can not however, account for the much larger surface ve¬ 
locity Vss ^ 4000 — TOOOmeVA reported in quantum os¬ 
cillations m- 

To confirm the results of our simplified analytical cal¬ 
culations, we have carried out a series of slab calculations 
with and without Kondo breakdown. As a first check, 
we simulated Kondo breakdown by setting the ampli¬ 
tude of the surface slave boson field to zero. Without 
surface Kondo breakdown (Eig 2(a)), the surface states 
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are heavy, but as expected, when the slave boson ampli¬ 
tude is suppressed to zero on the surface, the topological 
states become much lighter as shown in Fig. 2(b). More¬ 
over, the area of the light Fermi surface is significantly 
enhanced, corresponding to the extra density of carriers 
on the surface. 

Next, to confirm that the assumption of surface Kondo 
breakdown we performed a fully self-consistent mean field 
calculation in which both b{z) and A(^) were allowed to 
vary along the ^ direction. In this self-consistent calcu¬ 
lation, we chose parameters where the heavy and light 
bands cross at the three X-points, as seen SmBe mni]; 
from a technical stand-point, the crossing at the three 
X points is required to enhance the f-character, achiev¬ 
ing the Kondo limit Uf ^ 1 while still in the strong 
topological insulator phase. We achieve this crossing in 
our model by taking t = —f/2 , such that the heavy 
and light bands cross at the 3 different X points. We 
note that for the simpler choice t' = 0, bands cross 
at the F point, restricting the strong topological insu¬ 
lating behavior to the region to n/ <C 1. Fig 2 (c) 
presents the results of these calculations. We find as ex¬ 
pected, that the self-consistently determined slave boson 
amplitude b is depressed at the surface, corresponding 
to a reduction of the surface Kondo temperature. In 
our model calculation, for which T^ID ^ 0.025 found 
that b^ is depressed by 1/2 at the surface, corresponding 
to Tj^/Tk ~ 0.5 ^ {Tk/, a result consistent with 
equation (8). We also found that as the temperature is 
raised above the the surface Kondo temperature T^, the 
surface b collapses to zero, even though it is still finite in 
the bulk. The dispersion displays a corresponding tran¬ 
sition between heavy and light surface states as the tem¬ 
perature is raised through T^. In practice, the possible 
intervention of magnetism or other instabilities may per¬ 
manently prevent the low-temperature re-establishment 
of the Kondo effect at the surface. 

Chiral Kondo lattice: The consequences of surface Kondo 
breakdown are rather interesting. The decoupled f- 
electrons now interact with the chiral surface states to 
form a new kind of Kondo lattice in which the conduc¬ 
tion sea of chiral electrons is now singly degenerate, in¬ 
teracting with the local moments via a Hamiltonian 

Hckl = Y. ^c{k)clck + JY -Si + JaY^i- 

k i (ij) 

where J is the Kondo coupling, Sj is the decoupled local 
moment and where c\. creates a chiral, spin polarized 
surface states whose dispersion ec{k) is given in eq. 
The localized two-component electron field at site j. 



is constructed from the spin-polarized surface states by 
combining orbital and spin angular momentum. Jh ^ 


^t^f/cf is the antiferromagnetic Heisenberg exchange be¬ 
tween the local moments which is derived by a Schrieffer- 
Wolff transformation of eq. |S19[ The ground-state of 
the chiral Kondo lattice can be either magnetically or¬ 
dered state or a heavy fermion liquid where the local 
moments are screened by chiral fermions. As in conven¬ 
tional heavy fermion systems, near the quantum critical 
point that separates these two limits, there is the pos¬ 
sibility of strange metal behavior and superconducting 
ground-states. A 2D version of a similar model with ID 
edge states have been explored recently [34]. 

Even the magnetic phase diagram of this lattice is ex¬ 
pected to be rich. The RKKY interaction in a chiral 
Kondo lattice gives rise to Heisenberg, Dzyaloshniskii- 
Moriya and compass anisotropy terms [35| of strength J, 
D and A respectively. For AJ > 1 the ground state 
is an in-plane ferromagnet whereas for AJjD^ < 1, it is 
a spiral of variable pitch set by D/J|36]. Such a spiral 
state evolves into a skyrmion crystal in applied magnetic 
field [37]. Experimentally, in SmBe, indications of mag¬ 
netism have been observed as hysteresis in magnetoresis¬ 
tance experiments [38l l39] . The absence of ordered fer¬ 
romagnetic moment in XMCD experiments [25] suggests 
that the ground state might be a spiral with no net mo¬ 
ment. Magnetic order breaks the time reversal symme¬ 
try and can in principle gap the surface states. We note 
that of the most recent single crystals show no plateau 
in resistivity [25| indicating that magnetism may play a 
crucial role in the topological behavior. Better sample 
quality and tuning the chiral Kondo lattice through a 
quantum critical point opens up the possibility of exotic 
superconductivity on the surface of topological Kondo 
insulators. 

Thickness dependence of surface states: Perfect topolog¬ 
ical protection of the surface states requires an infinite 
bulk. When the thickness of the system is compara¬ 
ble to the decay length, f of the surface states, 

the two surface states on each surface mix, giving rise 
to a gap in the spectrum. Since the mixing is a zero- 
momentum transfer process, the gap opens at the dou¬ 
bly degenerate Dirac point as shown in Eig. 3(a). This 
effect has been observed in conventional topological insu¬ 
lators like Bi 2 Ses where the surface states get gapped for 
thickness 3-5 quintuple layers [40|. However with surface 
Kondo breakdown, the Dirac point shifts in the valence 
band, giving rise to a unique protection against decreas¬ 
ing thickness of the sample (Eig. 3(b)). Recent transport 
experiments on films about L=100-200nm thicknesses ob¬ 
serve topological surface states [41] . We predict that the 
topological states will persist even in ultra-thin samples, 
L ^ ^ lOnm, making StuSq a perfect candidate for 
thin film applications. This mechanism also gives rises to 
protection against small time reversal symmetry break¬ 
ing effects such as weak magnetism at the surface. 
Conclusion: To conclude, we have addressed the origin of 
the light surface states of topological Kondo insulators as 
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FIG. 3: (a) Without the surface Kondo breakdown, the sur¬ 
face states are susceptible to reducing thickness and and ef¬ 
fects of magnetism where an insulating gap starts developing 
around the Dirac point, (b) In the case of surface Kondo 
breakdown, the Dirac point sinks in the valence band, pro¬ 
viding an unusual protection to the surface states. The spec¬ 
trum stays metallic and qualitatively unchanged for moderate 
values of perturbation. 


a simple consequence of surface Kondo breakdown. Our 
results appear to reconcile the observation of light sur¬ 
face states and the peculiar shift of valence in Sm ob¬ 
served in X-ray measurements. This theory also predicts 
an increased stability of the the surface states against fi¬ 
nite thickness and magnetism may be important in thin 
film device fabrication. We have also argued that the 
interaction of the surface spins with spin-polarized light 
surface states is described by a chiral Kondo lattice: a 
two dimensional Kondo lattice with strong parity viola¬ 
tion, with the potential for a rich phase diagram in the 
vicinity of its magnetic quantum critical point. The de¬ 
tailed character of the surface states at low temperatures, 
with the possibility of a wide variety of magnetic or other 
orderings [42j [43] remains a fascinating topic for future 
research. 

Note added: We recently became aware of low 
temperature magneto thermo-electric transport 
measurements [44] which indicate the development 
of a very heavy surface Fermi liquid below 3K, with 
quasiparticle group velocities in the range 5 — lOmeVA, 
consistent with the development of a low-temperature 
surface Kondo effect. 
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Supplementary material for ^Kondo Breakdown In Topological Kondo Insulators’ 

This supplementary material presents the details of the surface state calculations. We first describe the Volovik 
method [ST] and then apply it to our model, with and without Kondo breakdown. 


VOLOVIK METHOD 

Surface states are evanescent solutions of the Hamiltonian which are nucleated by the boundary. The surface 
state wave function in the presence of a hard wall consists of an incoming and an outgoing wave moving in opposite 
directions whose superposition vanishes at the boundary as shown in Fig 1(a). Volovik’s methodfST] “unfolds” the 
reflection process into a transmission problem, with incoming and outgoing waves moving in the same direction 
(see Fig. lb). Similar “unfolding” techniques have been applied in wide variety of problems including the Kondo 
impurity model[S2]. The Volovik transformation is made by applying the reflection operator to the eigenvalue equation 
governing the incoming waves: under this a transformation, the odd parity terms in the Hamiltonian reverse sign at 
the boundary: in our particular, model the only odd parity term in the Hamiltonian is the hybridization, V{z), which 
in the Volovik mapping now changes sign at the boundary V{z) = Vsgn( 2 :) as shown in Figl(b). The Volovik method 
has the advantage that the ingoing and outgoing waves carry the same momentum; it also brings out the close analogy 
between surface states and Jackiw-Rebbi domain wall bound-states, with the additional technical advantage that it 
can handle cases where kz ^ 0. 




FIG. SI: Schematic of the bound state calculation (a) A hard wall with incoming and outgoing waves can be mapped onto (b) 
a reflection problem with a sign change in the hybridization. 

To be more specific, consider first a reflection boundary problem. The wavefunction consists of incoming and 
outgoing waves 


Reflection / 


= Elpi^2- 


(SI) 


Here V’i,2(-2) are eigenstates of the buik Hamiitonian defined on the haif-iine ^ > 0 and x is a spinor. The fuii 
wavefunction is given by 


ip{z) =ipi{z)-Ip 2 {z), {z>0). (S2) 

The boundary condition %p{z = 0) = 0 is satisfied by the fact that ^/’i(O) = ip 2 {Q)- 

The Voiovik method obtains an equivaient transmission probiem by appiying the parity operator to the eigenvaiue 
equation of the incoming wave PH+xl!i{z) = PExpi{z), to obtain 


H_'ipi{-z) = E'ipi{-z) 


(S3) 


where = PH+P^ is the parity reversed Hamiitonian. The corresponding transmission probiem is then given by 


Transmission 


V’ = -(pli-zWi-z) +'lp 2 {z) 0 {z) 
H = H_e{-z) + H+e{z) 


(S4) 
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Here 


(S5) 

is a symmetric domain-wall bound-state. In our particular problem the right and left hand Hamilonians are simply 
obtained by reversing the sign of the hybridization: H± = H{±V). 


CONTINUUM MODEL 


Our starting point is Hmf^ the mean field Hamiltonian, eq. 5 in the main text. We replace kz = {kp + idz) and 
expand Hmf to quadratic order in k so that 


H 


MF 


■ t{k^ 




1 0 

0 —a 


+ V {z){kx(Jx + + kz(Tz)Ti + A 


0 0 
0 1 


0(fc3), 


(S6) 


where 


V{z) = Usgn( 2 :). (S7) 

Here the r- index refers to the conduction/f-electron degree of freedom, so that = +1 describes a conduction 
electron while = — 1 describes an f-electron. The Pauli matrices (Tx^y^z carry the spin degree of freedom. We seek 
bound-state solutions of the form 




(S8) 


where k = {kx, ky^ kz). 

We will be interested in the small values of the transverse momenta, where \kx,y \ \kz\. This allows us to expand 

the Hamiltonian in in powers of k± Hmf = Hq Hi H 2 + where 

Ho = 2tkF(^^ d, + V{z)kFa,n+x(^^ +0{k^), (S9) 

Hi = V{z){kx(Tx + kyCFy)Ti, (SIO) 

H 2 = t{kl + kl) + 0{k^). (Sll) 

Our procedure is to first obtain the eigenvalues of the zeroth order Hamiltonian, and then to carry out 

perturbation theory on the remaining components Hi H 2 . This corresponds to first solving for the bound-state at 
kp = 0 . 

We need to treat positive and negative z separately. In the positive half-plane, {z > 0) the eigenvalue condition 
Hdip = gives 


E^ = 2tkF\^ dz^-{-VkpcrzTi^X ^ (S12) 

We may simultaneously diagonalize both energy and defining “up” and “down” eigenstates = Eip. For 

az = +1, the reduced eigenvalue equation is 


dzipt 


i { Ea VakF \ , 


(S13) 


where we have dropped the orbital index r for clarity. The solution for z > 0 can be written in the form of an 
evanescent wave. 


'Iptiz > 0) 


1 

AfVkF 


/l\ / iVkpOL 

V^/ cr “ xJaV‘^kp — 


^{ikz-H)z 


(S14) 
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where A/" is a normalization constant and kz , and r are defined as follows 


j r — aE 

~ 2takF ’ 

(S15) 

2r = E{\ T Of) — A, 

(S16) 

^/aV^kl - r 2 

2takF 

(S17) 


Similarly, the = — 1 solution is 


1 / 0 \ f - iVkFa 

^ NVkp V/cr y/aV‘^kp - 

Following the same procedure, the decaying solutions for z < 0 are 


^{ikz—K,)z 


'tp-f{z < 0 ) 

■(Pliz < 0 ) 


—iVkpOL 
Nvkp ^ + \/ 

iVkpd 


NVkp Vl/ o- + ^JaV‘^k\ — 


^{ikz+n)z 


^{ikz+K-)z 


(S18) 


(S19) 


Without Kondo breakdown 


We first consider the case the boundary conditions in the absence of Kondo breakdown. We need to match < 0) 
and > 0) at z = 0. This gives the condition: r = 0. We then find Eq^ k and kz for r = 0. 


p _ A 

(S 20 ) 

V 

K = - 

2t^/a 

(S 21 ) 

—Xry 

2t{la)akF 

(S 22 ) 


So far we have considered the // = 0 case. The solution for /i 7 ^ 0 can be obtained by 

H{X, li) = H{\ - M, 0) + /i. (S23) 


The remainder of the Hamiltonian is treated within perturbation 5H = Hi H 2 for small k^- The effective 
Hamiltonian derived by projecting 5H on to the surface states: = ('^al^^l'^cr'), where a, a' G are spin 

indices. The quadratic part of the Hamiltonian 5H cancels exactly: (' 0 cr|^ 2 |' 0 cr') = 0- The remaining linear part gives 

^Va 

Hu = = -2i^—{k,+iky), (S24) 

where the normalization = f dz{'ip^\xl>^) takes the form A/”^ = a{a + 1). The effective Hamiltonian is then 


= Ho+ 



ik+'\ 2 -y/aH 
0 ) a + 1 


(S25) 


where k± = k^ i iky and the corresponding dispersion is 


Eeff = Eo + 


a + 1 


(S26) 


where k± is defined as k']_ = k‘^ + ky. 
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Kondo breakdown 

In the presence of Kondo breakdown, the hybridization vanishes at the first site. As a result, c and / bands are 
subject to different boundary conditions. This effect can be captured by a simple phase shift such that the boundary 
condition > 0 )|^_q = < 0)|^_q is modified to 

(S27) 

Now the previous condition of r=0 is modified to r = ^/aVkp cos^. Hence the energy and the decay length takes the 
form 


K = 


K = 


X + 2y^VkF sin^ 


V 


2t^/a 


1 T 
sin 5. 


(528) 

(529) 


As before we consider the rest of Hamiltonian as perturbation 5H = Hi H 2 . The effective Hamiltonian is then 
H'^^, = ('0cr|^^|'0cr') However, the quadratic part {H 2 ) of 5H no longer cancels exactly and gives the diagonal part of 
the effective Hamiltonian 




where 


ak:^ l-exp(-2f^) ^ 


where 


= 


fc2 Q; + exp(-2f^) 

kZ + Ah 


leading to 


l-exp(-2f^) , 

= ta - 2 - 

a + exp-2f ^ 


While the linear part gives 


/ V ^.7 y ^ kz cos6 ^Ksin6 
Hfi = (V’tl77i|^-) = 2tk-—^/a 


kz'' 1+Q;exp(-2|^) 


The effective Hamiltonian is then 


77e// - Ho + 


HU Hn 


(530) 

(531) 

(532) 

(533) 

(534) 

(535) 


Combining all the pieces together we can diagonalise the effective Hamiltonian to get the energy: E^f f = Eq + ± 

|77ttl 


E, 


eff 


= E'n±2^^ 


kz 


A + aexp(-2f) 


-k I 


at' 


l-expi-2^) 
Q; + exp(-2f^) 


(S36) 
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